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— non-neutral  and  diamagnetic 
plasmas 
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Graduate  School  of  Frontier  Sciences,  The  University  of  Tokyo, 
Hongo,  Tokyo  113-0033,  Japan 


Abstract.  The  aim  of  this  paper  is  to  present  a unified  representation  of  relaxed 
states  generated  in  various  type  of  plasma  systems.  Self-organization  of  struc- 
ture is  the  central  them  of  the  recent  theory  of  relaxation  phenomena.  Invoking 
fluid  mechanical  models,  many  authors  have  discussed  various  relaxed  states  that 
constitute  a special  class  of  macroscopic  equilibria  (force-balance  states).  The 
present  theory,  however,  starts  from  a kinetic  model  of  relaxed  state,  which  en- 
compasses a variety  of  structures  including  both  electrostatic  and  electromagnetic 
configurations.  This  model  also  describes  the  macroscopic  relaxed  state  charac- 
terized by  the  Beltrami  and  (generalized)  Bernoulli  conditions,  which  gives  strong 
diamagnetism. 


I INTRODUCTION 

Plasmas,  viewed  as  a fluid,  can  produce  a variety  of  structures.  In  the  ideal 
(dissipation-less)  model  of  fluid  mechanics,  the  equilibrium  (force-balance) 
equations  may  have  infinite  number  of  solutions,  because  the  hyperbolic  part  of 
the  partial  differential  equation  demands  Cauchy  data  that  specify  the  profiles 
of  internal  fields.  To  select  a physically  realized  structure  from  the  total  set 
of  ideal  equilibrium  solutions,  we  invoke  the  effect  of  small  dissipations.  In 
approximately  collision-less  plasmas,  it  is  believed  that  the  effect  of  dissipation 
does  not  modify  the  equilibrium  structure,  while  it  chooses  a most  preferential 
profile  that  can  minimize  the  influence  of  the  dissipation.  This  stands  reason 
to  start  from  an  ideal  model,  and  then  search  for  a relaxed  state. 

Many  authors  have  discussed  various  relaxed  states  that  constitute  a spe- 
cial class  of  macroscopic  (fluid  mechanical)  equilibria.  The  pioneering  work  in 
this  context  was  J.B.  Taylor’s  model  of  the  RFP  (reversed-field  pinch)  plasma 
[1].  The  present  paper,  however,  presents  a kinetic  model  of  self-consistent  re- 
laxed states  that  are  characterized  by  some  rugged  constants  of  motion.  Both 
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electrostatic  and  magnetic  confinements  are  discussed  on  a unified  basis.  This 
model  also  derives  the  macroscopic  relaxed  state  characterized  by  the  Beltrami 
and  (generalized)  Bernoulli  conditions,  which  gives  strong  diamagnetism. 


II  KINETIC  EQUILIBRIA 

A Steady  states  in  kinetic  theory  and  relaxed  states 

We  begin  with  formulating  the  equilibrium  (steady  state)  in  the  phase  space 
of  kinetic  theory  that  is  spanned  by  x (coordinate)  and  p (canonical  momen- 
tum). Let  f(x,p,t)  be  the  distribution  function.  We  consider  a collision-less 
plasma.  The  evolution  of  / obeys  the  Liouville  equation  (Vlasov  equation) 

dj  + {//,/}  = 0,  (1) 

where  H is  the  Hamiltonian  of  a test  particle  moving  in  a mean  field,  and 
{ , } is  the  Poisson  bracket.  The  mean  field  A and  (j>  must  be  consistent  to  / 
through  Maxwell’s  equations.  We  consider  a stationary  state,  we  assume  that 
H is  independent  to  t.  Then,  a steady  state  for  (1)  satisfies 

{//,/}  = 0.  (2) 

The  stationary  equation  (2)  has  infinite  number  of  solutions.  Let  a;  be  a 
time-independent  quantity  that  commutes  with  //,  i.e.,  {H,  dj}  = 0.  The  H 
itself  satisfies  this  condition.  Suppose  that  we  know  N of  such  quantities  (con- 
stants of  motion).  Using  them,  we  generate  a distribution  (F  is  an  arbitrary 
smooth  function) 


/ = F(au-  ■ ■ ,aN).  (3) 

If  N is  equal  to  the  degree  of  freedom,  the  system  is  “integrable”,  and  (3) 
gives  the  complete  solution.  Our  motivation  in  this  work,  however,  is  to  find  a 
special  class  of  solutions  that  are  robust  (rugged)  against  various  perturbations 
destroying  the  invariance  of  fragile  constants  of  motion.  Only  a small  number 
of  them  are  robust  in  a sense  that  the  ensemble  averages  (or  the  total  sums) 
of  such  quantities  are  conserved. 

The  most  robust  steady  state  is  the  Boltzmann  distribution 

(4) 

where  Z (normalization  factor)  and  /?  (inverse  temperature,  or  a Lagrange 
multiplier)  are  positive  constants.  We  obtain  this  equilibrium  by  maximizing 
the  entropy  over  an  ensemble  that  is  characterized  by  a given  total  energy 
(i.e.,  a constant-energy  set).  If  vve  invoke  another  constant  G (in  an  ensemble 
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average  sense),  we  must  maximize  the  entropy  on  the  ensemble  defined  both 
by  H and  G.  With  restricting  the  totals  of  H and  G , we  obtain 

/ = Z~1e-0H-',G,  (5) 

where  7 is  the  second  Lagrange  multiplier.  Including  some  additional  con- 
stants of  motion,  we  obtain  a maximum  entropy  solution  that  is  more  re- 
stricted than  the  Boltzmann  equilibrium.  Such  a solution  remains  in  equilib- 
rium as  far  as  the  additional  constraints  work. 

B Equilibrium  with  momentum  conservation 

A symmetry  of  the  system  warrants  the  conservation  of  the  corresponding 
canonical  momentum.  Suppose  that  the  Lagrangean  L is  independent  of  a 
coordinate  a;0,  as  well  as  t.  Then,  the  Hamiltonian  H and  the  momentum 
p0  — dL/dx'0  are  conserved  ('  is  the  time  derivative).  With  an  arbitrary 
constant  C,  we  define 

H = H-  Cp0,  (6) 

and  consider  a distribution  [see  (5)] 

/ = Z~le~0**  = z~1e-0iH-Cp  o).  (7) 

This  solution  of  (2)  has  the  following  important  connotation  that  provides  a 
physical  meaning  for  the  parameter  C. 

When  we  discuss  a distribution  function  /,  we  consider  an  ensemble  of 
particles,  which  is  characterized  by  the  sum  of  the  Hamiltonian  over  the  all 
particles.  We  invoke  the  conservation  of  the  total  energy,  but  not  the  energy 
of  each  particle.  We  apply  the  same  framework  for  the  momentum  p0  in  (7). 
We  do  not  assume  the  conservation  of  H or  p0  for  each  particle,  while  we 
consider  an  ensemble  determined  by  the  totals  of  these  quantities.  Then,  the 
physical  meaning  of  H becomes  essential.  Indeed,  we  can  interpret  H as  the 
Hamiltonian  on  a moving  frame,  and  hence,  / = Z~le~0H  is  an  invariant 
of  the  collision  operator  (the  average  momentum  is  unchanged  by  collisions). 
This  robustness  of  po  warrants  the  use  of  p0  in  determining  the  ensemble. 

Let  us  revisit  the  change  of  variables  in  general  inhomogeneous  coordi- 
nate transform,  and  see  the  relation  between  H and  H.  Let  U be  a certain 
temporary-constant  velocity  field.  We  write  the  velocity  v of  the  laboratory 
frame  as 

v -v  + U,  (8) 

and  set  x'  — v.  The  Lagrangean  of  a charged  particle  (q:  charge,  m:  mass) 
can  be  written  as 


The  canonical  momentum  and  the  Hamiltonian  are 


P 

H 


dl 

— = m (i>  + U)  + qA  — mv  + qA, 
ov 


2 m 


qX\  - 


-V  +#  = -» 


rri 


u 2 + 


(9) 

(10) 


where 

~ 77? 

A = A + —U,  cf>  = 4>  - U ■ A.  (11) 

9 

The  effective  vector  potential  A includes  an  additional  term  that  yields  the 
Coriolis  force.  The  scalar  potential  <}>  has  received  the  (nonrelativistic)  Lorentz 
transform.  In  (10),  -mU2/ 2 is  the  centrifugal  potential.  The  transform  of 
the  Hamiltonian  and  the  momentum  can  be  written  as 


H = H -U  (mv  + qA)  = H -U  p,  (12) 

p = p = mv  + qA.  (13) 

Suppose  that  a component  p0  of  the  momentum  is  a constant  of  motion. 
Comparing  (6)  and  (12),  we  observe 

fl  = H — CPo  — H , (14) 

which  implies  that  H is  the  Hamiltonian  in  the  moving  coordinate,  and  C is 
the  velocity  (or  the  angular  velocity  if  we  define  p with  respect  to  an  angle)  of 
the  frame.  Hence,  the  solution  (7)  represents  the  “thermal  equilibrium”  in  the 
moving  frame.  Let  us  review  some  well-known  examples.  The  Hamiltonian  H 
includes  the  potentials  <j>  and  A that  must  be  self-consistently  determined  by 
the  field  equation  including  /. 

(1)  Electrostatic  equilibrium 

When  the  magnetic  field  produced  by  the  internal  current  is  negligibly  small 
in  comparison  with  the  externally  applied  magnetic  field,  we  can  use  the  elec- 
trostatic model;  A is  a prescribed  field,  while  <j>  is  determined  by  the  Poisson 
equation 


— A(j>  = -7i.  (15) 

Using  (10),  the  density  n is  given  by 

n(x)  = J f(x,v)dv  = noe~l3(~TnU*/2~qL’'A+q'1’) , (16) 

where  n0  = Z~l  J e-9mU/2 
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Here,  we  review  the  well-known  “thermal  equilibrium”  of  a cylindrical 
single-species  plasma  confined  in  a homogeneous  magnetic  field  ( B = Bez). 
Let  r-O-z  denote  the  cylindrical  coordinates.  By  the  symmetry  d/dO  = 0, 
the  canonical  angular  momentum  p$  = dL/dO'  = mrvg  + qrAg  is  conserved. 
Hence,  H — H — upg  (u>  = constant)  is  a constant  of  motion.  This  H is  the 
Hamiltonian  in  a rigid  rotation  frame.  Indeed,  setting  U = u>reg  is  the 
angular  velocity  of  the  rigid  rotation),  (10)  reads 

H = — v 2 - — (rw)2  - qrujAg  + q<f>  = H.  (17) 

The  equilibrium  f(H ) = Z~1e~PH  represents  a drift  Maxwellian  with  a con- 
stant angular  velocity  u.  Here,  we  seek  a solution  that  has  a constant  density 
inside  the  plasma,  i.e., 

k = vi‘  = b'p-‘,tf-  (18) 

The  vector  potential  for  the  homogeneous  longitudinal  magnetic  field  is  A — 
(rB/2)eg.  For  the  distribution  function  f(H)  = Z~le~l3ri  with  the  Hamilto- 
nian (18),  the  density  n is  constant  for  the  radius  r < a.  Then,  the  potential  is 
cj>  = — (qn/4c0)r2  [2].  To  satisfy  (18),  we  demand  — (m/2)(rw)2  — qrujAg  + q(f>  = 
0 [see(17)],  which  reads  as  the  familiar  equilibrium  condition 

U)2  + iOcU)  + -LO2  = 0,  (19) 

where  = qB/m  and  u2  = nq2/me0. 

(2)  Electromagnetic  equilibrium 

Next,  we  discuss  magnetic  confinement  solutions.  Let  us  consider  a slab 
plasma  with  B = Bz(x)ez.  We  assume  that  the  system  is  homogeneous  with 
respect  to  y ( d/dy  = 0),  so  that  the  canonical  momentum  py  = my'  + qAy(x) 
is  conserved.  We  consider  H = H — cpy  and  the  equilibrium  distribution 
f(H)  = Z~'e~PH , which  gives  the  drift  Maxwellian  with  a constant  velocity 
U = cey.  Here,  (10)  reads 

H - — v - —c  - qcAy  + q<t>=  H. 

We  note  that  this  H depends  on  x.  From  the  distribution  function  f{H), 
we  can  calculate  the  density  of  each  species.  For  appropriate  parameters, 
the  charge-neutrality  can  be  achieved,  and  hence,  (j>  can  be  set  zero  self- 
consistently.  The  drift  Maxwellian  yields  a finite  current.  The  magnetic  field 
produced  by  this  current  must  be  consistent  to  the  vector  potential  Ay(x)ey 
in  H.  Setting  Bz  — Botanh(x/f),  we  can  find  a self-consistent  solution  that  is 
called  the  Harris  sheet  [3].  The  magnetic  pressure  B 2 has  a dip  around  x = 0, 
where  the  pressure  is  confined.  This  equilibrium  is  based  on  the  conservation 
of  the  momentum  py.  Similar  calculations  in  cylindrical  coordinates  lead  to 
the  Bennet  pinch  equilibrium. 
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C Equilibrium  with  adiabatic  constants  and 
Beltrami-Bernoulli  conditions 

As  the  analysis  in  the  previous  section  shows,  a shear  flow  is  not  obtained 
from  the  model  of  momentum  conservation.  Here  we  examine  an  adiabatic 
invariant  model  to  find  an  equilibrium  with  a shear  flow.  Let  us  first  consider 
a certain  velocity  field  U and  define  a moving  coordinate  by  (8).  Here  we 
assume  that  U depends  only  on  x,  and 


U ■ Vx  = 0. 


(20) 


The  transform  of  coordinates  does  not  bring  about  temporal  variation  in  any 
field.  The  equation  of  motion  associated  with  the  Hamiltonian  H of  (10)  reads 


d__  _ _d_ 

dt  3 df>j 

_ _ d 
dt 


& = i fa  - > 


(21) 


= "4s  = ^ " A (~?u2 + #)  • (22) 


or,  in  the  form  of  Newton’s  equation  (with  dA/dl  — 0)  , 


d r~*  - . 1 w*  , 

m— v = q —S7<j>  fix(VxA)  + —S7 U 2. 


(23) 


In  (23),  the  Lorentz  force  includes  the  Coriolis  force.  We  define  a generalized 
magnetic  field  by  n = V x A.  We  assume  that  Li  is  a smooth  and  strong 
field. 

Here  we  introduce  an  essential  assumption 


fi  - aU, 


(24) 


where  a is  a certain  constant.  We  denote  _L  and  ||,  respectively,  the  perpen- 
dicular and  parallel  directions  with  respect  to  f 2.  For  “magnetized  particles”, 
the  generalized  magnetic  moment  jl  — Fx/I^l  's  an  adiabatic  invariant.  By 
(20)  and  (24),  we  have  Vy  = 0 for  every  scalar  field,  and  hence,  the  potential 
force  in  (23)  does  not  have  a parallel  component.  Therefore,  i>||  is  a constant 
of  motion,  and  gyration  average  of  the  kinetic  energy  mi2/ 2 does  not  change. 
Conservation  of  both  v2  and  fi||  yields  the  constancy  of  v\ , which,  together 
with  p.  =constant,  imply  that  |fi|  is  constant  (in  gyration  average)  through  the 
motion  of  particles.  By  (24),  ) t/|  is  also  constant  in  gyration  average.  Com- 
bining these  adiabatic  invariants,  we  observe  that  U ■ v = (v2  — f>2  — U2)  /2 
is  an  adiabatic  invariant. 

Using  this  adiabatic  constant,  we  define 

//  = //  — mU  ■ v = H — U ■ [p  - p),  (25) 
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where  p = mU  + qA.  We  may  rewrite  (25)  as 

mrr9  , m 
2^+*= 2 


» m_9  mTTn  . m.  ,rl9 

H = — u2  + — U2  Aq<t>=  — |w  - t/|2  + V5, 


where  <p  = rat/2/2  + qcj>  is  the  effective  potential.  The  equilibrium  define  by 
H reads 

f(H)  = = z~1e~^(m|”-Lr|2/2+v’),  (26) 


which  is  a drift  Maxwellian  with  a shear  flow  velocity  U.  We  note  that 
(24)  is  the  determining  equation  for  U,  which  corresponds  to  the  “Beltrami 
condition”  of  vortex  dynamics  systems  [4]. 

The  effective  potential  <p  determines  the  density; 

n = J f dv  = Hoe-11*  (n0  = Z~'  J e~^mi2/2  dv). 

The  pressure  p = nkT  ( kT  — 1//3)  satisfies  the  Bernoulli-like  condition 

Vp  + nV<p  = Vp  + nV(mf/2/ 2 + qcj))  = 0.  (27) 

We  note  that  the  Beltrami  and  Bernoulli  conditions  recover  the  fluid  mechan- 
ical equilibrium  model  of  diamagnetic  (high-beta)  structures  [4,5]. 


Ill  TOROIDAL  EQUILIBRIUM  OF 
NON-NEUTRAL  PLASMAS 

A Constant-density  equilibrium  in  a straight  annular 
geometry  with  sheared  magnetic  field 

In  this  section,  we  discuss  a more  specific  subject  aiming  at  development  of 
a new-type  of  non-neutral  plasma  trap.  We  consider  a toroidal  internal  ring 
system  that  can  produce  a magnetic  shear  configuration  [6].  To  understand  the 
basic  property  of  such  a configuration,  we  start  with  the  analysis  of  an  annular 
electron  plasma  column  (a  < r < b)  confined  by  both  longitudinal  ( Bz ) and 
poloidal  (Be  = BPa/r)  magnetic  fields  ( Bz  and  Bp  are  constants).  The  core 
region  of  r < a is  occupied  by  a conductor  that  carries  a longitudinal  current  to 
produce  the  poloidal  magnetic  field.  The  combination  of  Bz  and  Bg  produces  a 
magnetic  shear  configuration  with  spiral  field  lines.  The  vector  potential  of  the 
magnetic  field  is  A = ( (0,  ( Bz/2)r , —Bpa  log  r).  The  scalar  potential  <j>  is  given 
by  the  Poisson  equation  (15).  Here,  we  seek  for  a solution  with  a constant 
density.  Let  us  consider  a rigid  motion  of  the  plasma;  U — 1 (0,  ur,  Uz), 

where  u is  the  constant  angular  frequency  of  the  poloidal  rotation  and  Uz  is 
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FIGURE  1.  Toroidal  tlierinal  equilibrium  of  non-neutral  plasma. 


the  constant  longitudinal  velocity.  To  obtain  a constant-density  equilibrium, 
we  demand  -jU2  - qU  ■ A + q<p  = 0 [see  (10)],  which  reads 


—m 

~2~ 


q2nu 


2 , 

uj  u 4" 

m 2 mt0 


m . 


-f  q(BpaUz  + Ci)  log  r + qc2  - —U;  = 0. 


The  coefficients  of  r2  terms  yield  the  conventional  equilibrium  condition  (19), 
where  the  is  determined  by  the  longitudinal  magnetic  field  (u>c  = qBz/m). 
The  coefficients  of  logr  terms  yields  cx  — —BpaUz. 

For  this  equilibrium,  the  longitudinal  magnetic  field  Bz  plays  the  principal 
role  to  confine  the  plasma,  as  in  the  previous  cylindrical  equilibrium.  The 
poloidal  magnetic  field  Be  and  the  longitudinal  flow  Uz  produces  an  electric 
field  on  the  plasma,  while  it  does  not  work  to  confine  the  plasma.  Indeed, 
Bp  and  Uz  determines  the  coefficient  cx  of  the  electric  field,  but  they  are  not 
related  to  the  density  n0.  However,  the  poloidal  magnetic  field,  producing  a 
magnetic  shear,  has  an  important  function  to  stabilize  the  diocot.ron  (Kelvin- 
Helmholtz)  instability  [7]. 


B Generalized  thermal  equilibrium  in  a toroidal 
geometry  with  magnetic  shear 

In  a toroidal  system,  the  toroidal  symmetry  gives  the  constancy  of  the 
toroidal  canonical  angular  momentum.  We  denote  the  cylindrical  coordinate 
by  (r,  i?,z).  The  toroidal  angle  $ parallels  the  longitudinal  direction  z in 
the  previous  subsection.  The  poloidal  angle  is  no  longer  an  ignorable  coordi- 
nate, and  hence,  we  cannot  assume  the  conservation  of  the  poloidal  angular 
momentum  (which  played  the  essential  role  to  produce  the  constant  density 
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FIGURE  2.  Toroidal  equilibrium  of  non-neutral  plasma  base  on  the  drift  model. 


equilibria  in  Secs.  II B and  III  A).  Here  we  consider  a more  general  (non- 
constant  density)  solution  of  the  nonlinear  Poisson  equation  (15)-(16)  with 
setting  U — rfte#  (fl  is  the  angular  frequency  of  the  toroidal  rotation).  Writ- 
ing the  magnetic  field  as  B = Vt/>x  Vil  + B^rVH,  the  transformed  Hamiltonian 
is 

~ Til  Til 

H = H - = — v2  - y (rfi)2  - + q<t>- 

We  normalize  the  variables  as  x/A d = x,  and  (3qcp  — <j>,  where  A £>  is  the  Debye 
length.  The  nonlinear  Poisson  equation  now  reads 

— V2</>  — exp{— ^ + /3[m(rQ)2  / 2 + qQif>]}.  (28) 

When  n = 0,  (28)  is  just  the  Debye  shielding  equation.  The  plasma  flow 
(rotation)  has  the  effect  of  charge  neutralization.  To  confine  a high-density 
plasma  with  low  temperature  (Ad  -C  £),  O and  magnetic  field  must  be  carefully 
chosen,  because  e~ * becomes  a very  small  number,  and  the  density  tends 
to  localize  near  the  boundary,  producing  a sheath.  In  Fig.  1,  we  show  an 
example  of  thermal  equilibrium  with  Te  = lOeV  and  n0  = 10um  . When 
we  diminish  the  temperature,  equilibrium  density  becomes  strongly  localized 
near  the  boundaries. 

Another  model  of  non-neutral  plasma  equilibrium  considers  a drift  motion 
of  magnetized  particle.  For  a low  density  magnetized  non-neutral  plasma 
((Jc/wp  3>  1),  we  can  use  the  E X B approximation  of  the  flow  velocity  [8]; 
u = — V</>  x BIB2.  The  macroscopic  continuity  equation  in  a steady  state 
reads  V • ( nu ) = — [V</j  x \7(n/B2)\  ■ B — 0.  When  the  magnetic  field  has  a 
component  perpendicular  to  V0,  this  relation  demands  n = B2 F(<j>),  where  F 
is  a certain  smooth  function.  Using  this  n in  the  Poisson  equation  (15),  we  can 
find  equilibrium  solution.  In  comparison  with  (16),  this  model  can  produce  a 
larger  variety  of  solutions.  Figure  2 shows  typical  equilibrium  solutions  with 
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different  aspect  ratios  [9].  Because  the  drift  velocity  decreases  as  B increases, 
the  equilibrium  density  shifts  towards  a higher  field  region  (paramagnetism). 
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